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ABSTRACT 

Thermal conduction between a cool accretion disk and a hot inner corona can 
result in either evaporation of the disk or condensation of the hot corona. At 
low mass accretion rates, evaporation dominates and can completely remove the 
inner disk. At higher mass accretion rates, condensation becomes more efficient 
in the very inner regions, so that part of the mass accretes via a weak (initially 
formed) inner disk which is separated from the outer disk by a fully evaporated 
region at mid radii. At still higher mass accretion rates, condensation dominates 
everywhere, so there is a continuous cool disk extending to the innermost stable 
circular orbit. We extend these calculations by including the effect of irradia- 
tion by the hot corona on the disk structure. The flux which is not reflected 
is reprocessed in the disk, adding to the intrinsic thermal emission from grav- 
itational energy release. This increases the seed photons for Compton cooling 
of the hot corona, enhancing condensation of the hot flow and re-inforcing the 
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residual inner disk rather than evaporating it. Our calculations confirm that a 
residual inner disk can co-exist with a hard, coronally dominated, spectrum over 
a range of 0.006 < m < 0.016 (for a = 0.2). This provides an explanation for the 
weak thermal component seen recently in the low/hard state of black hole X-ray 
binary systems. 

Subject headings: accretion, accretion disks — black hole physics — X-rays: bi- 
naries — X-rays: stars 

1. Introduction 

It is well known that black hole X-ray binaries (BHXRBs) exhibit various spectral 
states. As recently reviewed by Remillard & McClintock (2006), such phenomena are quite 
common in BHXRBs. It is generally accepted that the transition from a spectrally soft to 
hard state is caused by a change in the accretion mode in the inner disk regions from a cool, 
geometrically thin structure at high accretion rates to a hot, geometrically thick, radiatively 
inefficient accretion flow (RIAF/ADAF) at low accretion rates. 

Many theoretical studies have attempted to identify the mechanism responsible for this 
transition, taking into account processes such as radial conductive energy transport (Honma 
1996; Manmoto & Kato 2000) and vertical evaporation (Meyer et al. 2000a; Rozanska 
& Czerny 2000; Spruit & Deufel 2002; Dullemond & Spruit 2005; Mayer & Pringle 2007; 
Bradley & Frank 2009). Of these, the disk evaporation model is one of the more promising. 
Thermal conduction between a cool disk and a hot corona can result in either evaporation of 
the cold disk or condensation of the hot corona (Meyer et al. 2000a,b; Liu et al. 2002, Liu et 
al. 2006, Meyer et al. 2007). At accretion rates higher than the maximal evaporation rate, 
fn > 0.03 (scaled by the Eddington rate M E dd = 1-39 x 1O 18 M/M g s" 1 ), for the specific 
case of a = 0.2, a full disk coexists with a corona. At very low accretion rates, fn < 0.006, 
the inner disk is completely evaporated into an ADAF. This gives an interpretive framework 
for features observed in BHRXBs, such as disk truncation as the cause of the spectral state 
transition (Liu et al. 1999; Meyer et al. 2000b; Liu & Meyer-Hofmeister 2001) and luminosity 
hysteresis between the hard-to-soft and soft-to-hard transitions (Meyer-Hofmeister et al. 
2005; Liu et al. 2005; Bradley & Frank 2009). 

It also predicts more complex behaviour at intermediate accretion rates, 0.006 < fn < 
0.03. Here the disk is truncated into a composite form with a coronal gap settling between an 
inner disk and an outer disk, with the inner disk carrying a fraction of the mass accretion rate, 
covered by a corona which carries the remaining material (see, e.g., fig.l in Meyer-Hofmeister, 
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Liu, Meyer 2009). This could explain the complex spectra seen in the intermediate state 
(Liu et al. 2006; Meyer et al. 2007), immediately following the soft state and predicts that 
a weak inner cool disk component can co-exist with a coronally dominated spectrum (Liu et 
al. 2007; Taam et al. 2008). 

This residual inner disk provides a potential explanation of the recent observation of 
a weak thermal component in some BHXRBs in the low/hard state (Miller et al. 2006a,b; 
Tomsick et al. 2008, Reis et al 2009; 2010; Chiang et al 2010). However, the calculations 
to date have only studied evaporation as a means of coupling the disk and corona, whereas 
radiative coupling from irradiation of the disk by the corona is also important. The corona 
illuminates the disk, and most of this irradiation is reprocessed, enhancing the luminosity 
of the disk, and these photons provide additional Compton cooling of the corona so there is 
strong feedback between the disk and corona (Haardt & Maraschi 1993). Here, we investigate 
the effect of this additional radiative coupling between the disk and corona, especially to 
determine whether the inner disk can survive or is completely evaporated by irradiation. 

In §2 we briefly introduce the disk corona evaporation/condensation model and its 
extension to incorporate the effect of coronal irradiation. Numerical results from these 
models are presented in §3, showing the strength of inner disk relative to the corona/ADAF. 
Our conclusion is presented in §4. 

2. The model for an inner disk around a black hole including irradiation 

The disk/corona model adopted here is based on the studies of Liu et al. (2006) and 
Meyer et al. (2007) with modifications recently incorporated in the works of Liu et al. 
(2007) and Taam et al. (2008). It is assumed that an ADAF-like corona (given by the 
self-similar equations of Narayan & Yi 1995, with a = 0.2 and = 0.8) lies above a thin 
disk. The corona is heated by the viscous release of gravitational energy of accreted gas and 
dominantly cooled by vertical conduction and Compton scattering of soft photons emitted 
by the underlying disk. In the vertical transition layer between the disk and corona, an 
equilibrium is established between the conductive flux from the upper corona, bremsstrahlung 
radiation, and vertical enthalpy flux. For a given distance from the black hole, a fraction of 
the disk gas is heated and evaporated to the corona when the conduction flux is too large 
to be radiated away. On the other hand, a certain amount of coronal gas is cooled down, 
condensing to the disk if the bremsstrahlung radiation is more efficient than the conduction. 
At accretion rates around a few percent of Eddington value, gas evaporates from the disk 
to the corona around a few hundred Schwarzschild radii and partially condenses back to the 
disk in the innermost region (Liu et al. 2006; 2007; Meyer et al. 2007). 
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2.1. The condensation feature without irradiation 



In the case of conduction dominant cooling the condensation rate is given as (Liu et al. 
2006; 2007; Meyer et al. 2007; Taam et al. 2008) 

m cnd (r) = 3.23 x 10~ 3 a~ 7 m 3 [l - 6(r/r d ) 1/2 + 5(r/r d ) 3/5 ] , (1) 

where r d is the critical radius from which gas flowing changes from evaporation in the outer 
region to condensation in the inner region. Here, r d , is expressed as 



r d = 0.815o~ 28/3 m 8/3 . 
In the case of Compton dominant cooling the condensation rate is 
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and r d is the critical condensation radius, 
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T e fj max is the maximal disk temperature, which depends on the condensation rate integrated 
from outer edge of the inner disk (r d ) to the distance r tmax = 49/12 



T eS>max = 1.3348 x 10 7 m- 1 / 4 m;/ 4 (r tmax )K. 
The radial distribution of disk temperature is approximately expressed as, 
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For a black hole with given mass m (in unit of solar mass), accretion rate rh and 
viscosity parameter a, the condensation rate can be determined by EqQ] or Eqj3l In the 
Compton dominant case, an effective temperature (T eff max ) is presumed for calculating the 
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condensation rate, with which a new temperature is derived from Eqd Iterative calculations 
are carried out until the presumed temperature is consistent with the derived value. Compton 
cooling becomes dominant when the accretion rate (m) is high or/and the disk is heated to 
a high temperature (for details see Liu et al. 2007). At low accretion rates, Compton 
cooling is either negligible or only dominant in a very narrow inner region. In this case, the 
total condensation rate (and hence the disk temperature) is an integral over the conduction 
dominant region and Compton dominant region (if exist), and the size of the inner disk is 
determined by Eqf2J 



2.2. Inclusion of irradiation 



Since the inner disk is optically thick, the irradiating photons are primarily reprocessed 
as blackbody radiation. These additional soft photons propagate through the corona and lead 
to more efficient Compton scatterings. This leads to enhanced coronal cooling, so a greater 
amount of coronal gas condenses into the inner disk. Therefore, the effect of the irradiation 
is to increase the condensation, building up a relatively strong inner disk. This increase is 
dependent on the irradiative luminosity, the absorption/scattering of the associated photons 
by the inner disk, and the covering factor of the inner disk as seen from the ADAF. The 
value of the albedo is usually taken as a constant equal to 0.15 (Haardt & Maraschi 1993). 
Both the irradiation luminosity and covering factor are determined by the accretion rate, 
the latter dependence resulting from the dependence of the inner disk size on the accretion 
rate. 



2.2.1. The covering factor of a disk to an ADAF 



Assuming that the radiation from the ADAF radiation can be represented as a point 
source lying above the disk center at a height of H s , the covering factor of a disk ring at 
distances between R in < R < R ut is given by 
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For a disk characterized by R lTL = and R out = oo, Eqj9] yields / = 1/2, indicating that 
half of the ADAF radiation is intercepted by a full disk. For a small disk ring characterized 
by 3R S < R < 10R S , the covering factor is 0.2, provided that the height of the irradiation 
source is, H s = i? in = 3R S (R$ is the Schwarzschild radius). Therefore, only a small fraction 
of the ADAF radiation illuminates the truncated inner disk ring. 
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2.2.2. The irradiations from the corona 



Assuming that the total intrinsic luminosity of the corona above and below the disk is 
given as L C) i n , the illumination flux to the disk surface per unit area is, 
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where a fraction of the irradiation flux, a, is assumed to be reflected at the disk surface. To 
simplify the calculation of the condensation/evaporation rate, Eq{TU]is approximated as, 
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which is equivalent to assuming the covering factor per unit surface area, 4^2+^2) 
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The coefficient of 3 in Eq.l flTT) is a normalization factor determined 
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by the fact that a full disk covers half of the sky of a corona as shown in §2.2.11 (that is, 
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to the exact expression for the irradiating flux given by eq.l fTOl) at R = (49/36)i? s where a 
standard disk reaches its maximal temperature and, hence, is a good approximation. 

The illumination flux heats the optically thick disk and is eventually re-emitted as soft 
photons. Thus, the radiative flux from the disk is composed of both accretion released energy 
and irradiative energy by the corona, which we express as 
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The corresponding disk temperature is modified by irradiation as, 
T eff (r) = 2.05T eff , 
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Therefore, the disk temperature is raised by a factor of 
added to the expression of T c fj imax (EqJTj) in calculating the condensation rate with EqsjSJEl 
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This factor is 



3. Results 



Given the mass of a black hole and an accretion rate, the condensation rate, corona 
luminosity and size of an inner disk can be calculated. From these quantities, the covering 
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factor and, hence, the irradiating flux at any given distance can also be determined. As 
the density of the soft photons, representing reprocessed irradiating photons and photons 
originating from the viscous dissipation in the disk, differs from the unirradiated case, iter- 
ative calculations are performed until the derived soft photon density is consistent with the 
presumed one. For M = 1OM , a = 0.2, H s = 3R S , a = 0.15, m = M/M Edd < 0.03 (so 
that a composite corona/ADAF geometry can form) we calculate the condensation of this 
coronal flow to an inner, residual disk and determine the relative strength of the corona with 
respect to the disk. 

Fig. [1] shows the electron temperature in the corona as a function of mass accretion 
rate. There is a marked change in behaviour at m = 0.015, with irradiation making very 
little difference below this point. This is because at low mass accretion rates the cooling in 
the corona is dominated by vertical heat conduction rather than by Comptonisation. Thus 
increasing the soft photons from the reprocessed radiation flux makes little difference to 
the coronal temperature, so the condensation rate is unaffected. However, at higher mass 
accretion rates, Compton cooling of the corona becomes important and the inner disk is 
large. The inner disk subtends a large angle to the corona radiations, further increasing 
irradiation. Thus the seed photon luminosity largely increases, leading to a decrease in 
electron temperature of the corona. This decreases the conductive heating to the transition 
layer, and so leads to an increase of the condensation rate. We show this in Fig. defined as 
the integrated condensation rate from the outer edge of the inner disk to the innermost stable 
circular orbit, ISCO, (so the mass flow rate through the innermost part of the corona at the 
ISCO is m — rh cn d). However, this increased condensation has very little effect on the size of 
the residual inner disk (Fig. [3]) as by this stage the outer edge of the inner disk is already 
fairly large. Irradiation is preferentially of the inner disk so the increased condensation adds 
material to the inner disk rather than increasing its outer radius. 

While irradiation of the disk makes little difference to the corona at the lowest mass 
accretion rates at which the composite geometry can exist, it does make a difference to 
the observed emission from the inner disk because the intrinsic accretion-heated inner disk 
emission is extremely dim compared to the corona illumination. FigJH shows the inner disk 
temperature for the irradiated and non- irradiated case. Irradiation heating means that the 
disk has a higher temperature (and luminosity) than expected from its mass accretion rate 
alone. This can be compared to the compilation of observations of the residual thermal emis- 
sion in BHXRBs of Reis et al (2010), where the typical temperature of the soft component 
is ~ 0.2 keV. There is probably a colour temperature correction of ~ 1.6 — 1.8 between the 
observed and effective temperature (Shimura & Takahara 1995), so the data probably corre- 
spond to an effective temperature of ~ 0. 10 — 0.15 keV, matching well with the predictions 
of the irradiated inner disk for low mass accretion rates (rh < 0.02). 
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The enhanced luminosity of the disk also impacts on the observed spectrum of the 
corona. At the lowest luminosity the coronal radiation is dominated by bremstrahlung 
rather than Compton scattering, so the (small) increase in seed photon flux makes little 
difference to the spectrum. However, Comptonisation increases in importance as the mass 
accretion rate increases, and is comparable to bremstrahlung from the transition layer at 
rh = 0.015. Above this, Comptonisation is dominant, and the ratio of (reprocessed) disk to 
coronal emission sets the slope of the Comptonised spectrum so that the larger this ratio, 
the softer the spectrum (Haardt & Maraschi 1993). 

Most of the spectra in Reis et al (2010) show a weak soft component co-existing with 
hard spectra, with T < 1.8. From Compton cooling models, both numerical (EQPAIR: Coppi 
1991) and analytic (Pietrini & Krolik 1995) require that the ratio Ld/L c < 0.3. This ratio 
can be estimated analytically as a function of the covering fraction of the disk, /, and its 
albedo, a, so that L&j L c = f(l — a)/ [1 — /(l — a)} (Haardt & Maraschi 1993). This is a 
lower limit as there should be some intrinsic dissipation from the material accreting through 
the disk. Hence T < 1.8 (equivalent to 0.3 > Ld/L c > /(l — a)/ [1 — /(l — a)]) requires 
/ < 0.27 for a = 0.15. Such a low covering fraction requires that the residual inner disk 
is rather small, corresponding to R out = 13.5R S . Such a small disk is only possible at very 
low mass accretion rates, so there is very little range in rh over which the residual disk can 
co-exist with a hard spectrum. 

However, this conclusion is very dependent on the assumed albedo. A more reflective 
disk (so that less of the irradiating luminosity is reprocessed) can give harder spectra for 
larger inner disks. However, hard spectra peak at high energies, where Compton down- 
scattering means that the energy cannot be completely reflected. The leads to a maximum 
possible albedo of a ~ 0.6, corresponding to a covering fraction / < 0.58. This estimate 
appears not to constrain the disk size. We note that the disk accretion could contribute a 
luminosity comparable with the reprocessed disk luminosity in the case of a large albedo, 
a = 0.6. Including this effect, the covering factor required to produce spectra T < 1.8 is 
smaller. 

We confirm these estimates with the numerical results of L&f L c from our models. This 
ratio is shown in Fig. [5] for the representative albedos of 0.15 and 0.6, and the requirement 
that the ratio is less than 0.3 in order to produce hard spectra (r < 1.8) directly constrains 
rh < 0.01 or 0.016 for a = 0.15 and 0.6, respectively. 

Thus the model can indeed produce a residual inner disk temperature and hard coronal 
luminosity to match the observations of the small thermal component seen together with a 
hard X-ray continuum, but only for a limited range in mass accretion rate of 0.006 < rh < 
0.01 for mostly neutral reflection or 0.006 < rh < 0.016 for a highly ionised disk. Thus the 
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model predicts that the composite geometry (outer disk, coronal gap, with mass accretion 
rate in the innermost regions split between a weak disk and coronal flow) can only be seen over 
a range in luminosity of at most factor ~ 7 assuming an efficiency of the coronal flow, r] oc rh. 
Even this maximum range in luminosity is slightly smaller compared to the observations of 
Reis et al (2010). These show the thermal component together with hard spectral index 
(r < 1.8, so J1650-500 with T = 2.1 excluded) over a range of about a factor 10 in Lx/L^dd 
(more if XTE J1118+480 is included), but we note there are substantial uncertainties in 
system parameters, so this can be consistent with the data. A bigger discrepancy is that 
the mass accretion rate of rh = 0.006 — 0.016 in the model corresponds to a luminosity of 
L/LEdd ~ 1 — 9 x 10~ 3 (assuming a self-similar ADAF with energy conversion efficiency 
rj — 0.1 x jj^_) , somewhat smaller than that seen in the data which cluster around L/LEdd ^ 
5 — 9 x 10~ 3 (Reis et al 2010). However, the scaling from mass accretion rate to luminosity 
is dependent on the model parameters chosen, particularly on a (Esin et al 1998). 

At higher mass accretion rates, the spectrum softens strongly as the increased seed 
photon flux leads to increased Compton cooling of the corona. This means that the spectral 
transition is triggered earlier than expected from the statement that the disk should exist 
at all radii. Models with rh ~ 0.02 do not have a continuous disk. There is still a coronal 
gap, but the spectrum is formed mainly in the innermost regions, and at these radii the disk 
covers almost half the sky as seen from the corona. Comptonised emission calculated from 
such a composite accretion flow with electron temperature and optical depth determined by 
our model for rh ~ 0.02 yields a photon index T ~ 2. Thus the transition from the low/hard 
to intermediate state occurs at rh ~ 0.02 or L/LEdd ~ 1-3 x 10 -2 . 



4. Conclusions 

We investigate the influence of coronal radiation on the properties of an existing residual 
inner disk. For accretion rates close to the hard-to-soft transition rate (m ~ 0.03), the disk 
is only truncated in a very narrow radial extent by mass evaporation. In this case, the disk 
is described by a geometrically thin inner region, a narrow geometrically thick intermediate 
region, and a geometrically thin outer region. Here we find that including irradiation results 
in efficient Compton cooling and condensation of coronal gas to the inner disk, with the 
emergent spectrum much softer than predicted without reprocessing. This means that the 
transition from the intermediate to low/hard state occurs at a lower mass accretion rate, at 
rh ~ 0.016 — 0.02 for the specific parameters assumed here, rather than at rh ~ 0.03 without 
irradiation. 

Below this, evaporation from the outer region results in an increased radial extent of the 
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ADAF region and a smaller inner disk. The resulting geometry leads to a greater fraction 
of photons which can escape from the ADAF without irradiating the inner disk. Therefore, 
at sufficiently low accretion rates (m < 0.015), the irradiation of the inner disk is small and 
the resulting Compton cooling does not significantly influence the dynamical process in the 
corona. Consequently, the condensation rate is not affected, while the inner disk is heated 
to a higher effective temperature by irradiation. 

Thus a weak, inner disk can exist in the low/hard state even if irradiation is taken 
into account. The ratio of luminosities emergent from the disk and corona decrease with 
decreasing mass accretion rates, until fn ~ 0.006 at which mass condensation no longer 
occurs and the inner disk cannot be sustained. Our results are in broad agreement with 
recent X-ray observational results that systems with a weak disk can exist in the low/hard 
spectral state (Reis et al. 2010). 
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Fig. 1. — The variation of the electron temperature in the corona with respect to the mass 
accretion rate for the cases with irradiation (solid line) and without irradiation (dotted 
line). Irradiation does not affect the temperature at low accretion rates, while it leads to 
lower temperatures at higher accretion rates. The maximum temperature corresponds to 
the transition where the dominant cooling process changes from conduction at low accretion 
rates to Compton scattering at high accretion rates. 
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Fig. 2. — Comparison of the condensation fraction of the coronal accretion flow for the cases 
with irradiation (solid line) and without irradiation (dotted line). The condensation fraction 
is unaffected by the irradiation at low accretion rates, while it is increased by the irradiation 
at high accretion rates. 
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Fig. 3. — The size of the inner disk is shown as a function of accretion rate. The two curves 
including irradiation and excluding irradiation overlap, indicating that there is no influence 
from irradiation. 
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Fig. 4. — The maximal effective temperature of the inner disk heated by accretion and 
irradiation is shown as a function of the mass accretion rate. The solid line corresponds to 
the case with irradiation and the dotted line for the case without irradiation. 
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Fig. 5. — The ratio of the luminosity from the disk relative to the corona is shown as a func- 
tion of the mass accretion rate. The solid line corresponds to the case with irradiation and 
the dotted line to the case without irradiation. Here L c denotes the emergent coronal lumi- 
nosity where the fraction absorbed by the inner disk is deducted from the coronal emission. 
The luminosity ratio is significantly increased by the irradiation, however, it steeply drops 
to zero when the accretion rate decreases to the critical value of ~ 0.006. The influence of 
the albedo is shown by the dashed line, where a = 0.6 is assumed. The effect of irradiation 
is much reduced if a larger fraction of the irradiation flux is reflected at the disk surface. 



